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1 Introduction 



Recently, there has been considerable interest in non-hermitean random matrix theory 
||T| as applied to a number of interesting physical situations. A particularly interest- 
ing issue is that of the localization-delocalization transition in non-hermitean random 
Hamiltonians. The earliest work to address this issue appears to be that of Miller 
and Wang in their study of a problem in fluid flow. More recently, Hatano and 
Nelson have mapped the problem of the vortex line pinning in superconductors to 
a problem involving a non-hermitean random Hamiltonian. The prototypical Hamil- 
tonian 

H = Ho + W (1) 
is the sum of a deterministic non-hermitean hopping term 

^oii = ^ (e^ 5m,i + e-'' , ^,J = 1,■■■A^ (2) 

(with the obvious periodic identification i + N = i of site indices) describing the 
drift of a vortex line being driven by an external current and the hermitean random 
potential term 

Wij = Wi 6ij . (3) 

describing the pinning of the vortex line by impurities. Here the real numbers Wi 
are drawn independently from some probability distribution P (w) (which we will 
henceforth take to be even for simplicity.) The number of sites N is understood to 
be tending to infinity. We will also tak t and h to be positive for definiteness. This 
problem has been studied by a number of authors [H, H, H, 0, |^ following Hatano 
and Nelson. Note that the Hamiltonian not only breaks parity as expected but is 
non-hermitean, and thus has complex eigenvalues. It is represented by a real non- 
symmetric matrix, with the reality implying that if E is an eigenvalue, then E* is 
also an eigenvalue. In this paper, we extend and generalize the analysis in [^, which 
we will refer to extensively as FZ. We will review some of the salient points in FZ 
briefly; for further details, the reader is referred to FZ and to the literature. 
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With no impurities {wi = 0) the Hamiltonian is immediately solvable by Bloch's 
theorem with the eigenvalues 

E„ = t cos -ih), (n = 0, 1, ■ ■ ■ , iV - 1) , (4) 

tracing out an ellipse. The corresponding wave functions ■?/']"'■' ~ exp 2irinj/N are 
obviously extended. 

With impurities present, "wings" emerge out of the two ends of the ellipse. (See, 



for example, figure 1 of FZ. See also (^, 0) below.) Some eigenvalues are now real. 



Evidently, the "forks" where the two wings emerge out of the ellipse represent a non- 
perturbative effect, and cannot be obtained by treating the impurities perturbatively. 

A simple argument shows that all localized eigenstates of (|l]) correspond to the real 
eigenvalues on the wings. Conversely, the states corresponding to complex eigenvalues 
are extended, that is, delocalized. In this sense, non-hermitean localization theory is, 
remarkably enough, simpler than hermitean localization theory. 

One way to think about this problem is to imagine starting with Hq and then 
increasing the randomness (measured by a parameter 7, say) and ask if there is a 
critical strength of randomness 7ci at which localized states first appear. There may 
also be another critical strength 7c at which all the extended states become localized. 
Obviously another way of thinking about the problem is to start with the hermitean 
random problem {h = 0) in which according to Anderson and collaborators all 
states are localized. We then increase the non-hermiticity and ask about the critical 
strength of the non-hermiticity at which extended states first appear. Similarly, 
we can ask, for given h and 7, the critical hopping amplitude tc needed for the states 
to become de-locahzed. 

In FZ, the emphasis is to understand analytically the essential physics involved, 
including the non-perturbative emergence of the two "wings" along the real axis. 
While some general analytic statements can be made, exact results are obtained in 
FZ only in two simplifying limits: the single impurity limit and the one way limit. 

In the single impurity limit, the A^ random impurities in (^) is replaced by a single 
impurity; in other words. 



Wij = w6i,i6j^i (5) 

It was argued in FZ that this simphfication, though drastic, captures some of the 
essential physics. 

In the one-way limit, the parameters in (^) are allowed to tend to the (maximally 
non-hermitean) limit h ^ oo and t — > such that 

te'' ^2 (6) 

The spectrum (^) changes into 

E„ = exp^, (n = 0,l,---,iV-l), (7) 

and the ellipse associated with expands into the unit circle. 
As is standard, we are to study the Green's function 



G{z) 



z-H 

The averaged density of eigenvalues, defined by 



p(x, ^ (1 ^ 5{x - ReE,) 5{y - ImE,)) , (9) 



is then obtained as 



with z = X + iy. Here we define 



1 



G{z)^{-tTG{z)) (11) 

and < . . . > denotes, as usual, averaging with respect to the probability ensemble 
from which H is drawn. A small subtlety is that while G{z) is ostensibly a function 
of z only, {G{z)) depends on both z and z*. A careful discussion is given in 0. 
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It is also useful to define the Green's function Go{z) = jzj]^ the absence of 
impurities and the corresponding Go{z) = ti jzjj^- For the spectrum in (H), we 
can use Cauchy's integration to evaluate 



Go{z) 



de 1 



271 z — t cos(6 — ih) 
2 r dv 1 



t J 2m {v—v+){v — v. 



(12) 



where we have changed variable to t; = e^^^^ and v± = j ± J ^ — 1. We obtain 



1 



Go{z) = -^^^ + I -te') (13) 

(There is ostensibly another term involving 9(\ z — \J z^ —t'^ \ — t(^\ To show that 
this term can be dropped and that the boundary defined in the step function in ( [T3| ) is 
just the ellipse in (^), write z = t cos(mi — ^2) = t(cos Ui cosh U2 +i sin Ui sinh U2) 
with < Ml < 27T and U2 > 0. To cover the complex plane, we can restrict U2 to 



be non-negative. We define a/z^ — = t(cos ui sinh U2 +i sin ui cosh U2) so that 



when z is in the upper half plane V-z^ — is also in the upper half plane. With 
these conventions, the second step function can be dropped. Then the step function 
in ( p!3D implies U2 > h. At U2 = h, z = x + iy traces out the ellipse 



+ — 

(cosh h)'^ (sinh hy 



+ T^^^ = (14) 



in agreement with (^. Note that this goes to the appropriate function in the her- 
mitean limit h — > 0. In the one way limit, (|13]) simplifies to 

1 



Go{z) = -ei\z\-l) (15) 
In FZ, the following probability distributions were studied: the sign distribution 

P(w) = ^[5{w-r) + 5{w + r)] (16) 
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with some scale r, the box distribution 



(17) 



obtained by "smearing" the sign distribution, and the Cauchy distribution 



P{w) 



7 



1 



(18) 



TT W 



,2 _^ y 



,2 ' 



with its long tails extending to infinity. 

We might also wish to consider the effect of diluting the impurities by setting 
randomly some fraction d of the w^'s to zero. In other words, given P{w) we can 
consider 



with < d < 1. As d increases, there should be more extended states. 

The problem of determining the spectrum of H with the Cauchy distribution was 
solved explicitly in FZ in the one-way limit. The density of eigenvalues is made of 
two arcs and two wings. (For completeness, we will derive a generalization of this 
result later in this paper.) The complex eigenvalues form the union of two arcs of a 
circle of radius one. More precisely, for y 0, 



The upper arc is that part of a semicircular arc of a unit circle (centered at the 
origin) that remains in the upper half plane after being pushed a distance 7 downward 
along the imaginary axis. (The lower arc is of course the mirror image of the upper 
arc.) Each of these arcs is thus of length arc cos (27^ — 1) < tt and carries narc = 
{N/27r) arc cos (27^ — 1) < (A^/2) eigenvalues. This means that the arcs exist only 
as long as 7 < 7c = 1. 

The rest of the eigenvalues, which do not have space to live on the arcs, must 
have snapped onto the real axis and formed "wings". The eigenvalue density along 
the real axis is given by 




(19) 




(20) 
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Pwing(a;, y) = ^ ^— ^ Oix" + 7' - 1) 5{y) (21) 

We thus see that for 7 < 7c = 1, there are two wings that bifurcate from the arcs at 
X = ±Xc = ±-\/l — 7^. Integrating over ( PT| ) we find that the fraction of eigenvalues 
that reside in the wings is (n^ings/N) = 1 — (2/7r) arc cos 7 = 1 — (l/vr) arc cos (27^—!), 
which together with the fraction of eigenvalues that reside in the arcs, sum up exactly 
to 1. As 7 tends to 7c, Xc becomes smaller, and vanishes at 7 = 7c. At this point 
the two wings touch at the origin and the two arcs disappear completely. All the 
eigenvalues have snapped onto the real axis. 

Note that in this problem the critical coupling 7ci defined earlier is equal to 0^. 
Localized states appear no matter how small 7 may be as long as it is non- vanishing. 
(Since we have already gone to the maximally non-hermitean limit we cannot ask 
about he-) 
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2 Multiple Scattering and "Two Way" Hopping 



In this section we go beyond FZ and discuss what we can obtain analytically without 
taking the "one way" limit in (^. With the particle now able to hop both ways, the 
calculation of G{z), which as usual can be expanded formally as a series: 

G{z) = Go{z) + Go{z)WGo{z) + ...,, (22) 

now involves the combinatorial problem of summing over all possible histories in 
which the particle, starting from some intitial site i, propagates to the site j, scatters 
off the impurity there, then propagates to the site k, scatters off the impurity there, 
and so on, until it arrives at some final site /. This is precisely what is expressed by 
M) which we can rewrite as 



G 



, 1 - E WkGoPk I 

\ k I 



Go (23) 



Here we have defined the projection operator 

Pk = \k){k\ (24) 

In this section we will not commit to any specific form of Hq. We will merely 
assume that Hq is translation invariant. One can average (^) with one's favorite 
P{w). In practice, however, it is more useful to carry out some partial summation 
over the effects of repeated scattering on a specific impurity. Let us define 



Vk = 



1 - Wk{k I Go I A;) 
Wk + wl{k I Go I A;) + wl{k \Go\kf +■■■ (25) 



and 



= VkGoPkGo (26) 
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By translation invariance, {k \ Gq \ k) = Gq{z). The matrix {gk)ij describes the 
histories consisting of propagation from some initial site i to the site k, scattering off 
the impurity at site k, once, twice, . . ., any number of times, and finally propagating 
to some final site j. In other words, only scattering off the impurity at site k is 
allowed. Next we define 

Gk=gk + 9kGo'T.Gi (27) 
Then the desired Green's function is given by 



G — (jQ + Cfc 



(28) 



Putting these equations together, we obtain 



G 



VkGoPk 



n -1 



Gn 



(29) 



It is straightforward to expand this expression and show that it leads to (p3D, of course. 
This expression, while equivalent to (PBj), is however more suitable for averaging. 

For any arbitrary P{w) it is presumably not possible to average (^) explicitly 
and obtain {G{z)) in closed form, not any more than it is possible to obtain {G{z)) 



in closed form for the hermitean problem. Indeed, (|2^) was derived by only assuming 
translation invariance for Hq, whether hermitean or non- hermitean. 

It is, however, always possible to expand ( P^ ) to any desired powers of v, and 
average. To write down the result, it is best to define Dyson's irreducible self energy 
Y.{z) by the standard formula 



To third order in v we obtain 



(30) 



E(^) 



(v) 



{v)Go 



Gl 



— 2 

Gl 



(31) 
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Here we have defined 



and 

= {k\Go\k) = ^trGo (33) 

and 

Gl = {k I Gl I A:) = ItrGg = - (34) 

For Ho in (|I]), we can read off Go and lience Gg from (|13|). Wliile 1^(2:) is in general a 
matrix, it is merely a number to this order. Thus, to this order we obtain a remarkably 
simple result 

Gi^ = Goiz-Eiz)) (35) 

with J2{z) given in (|3TD . Thus, the reader can simply compute (v) with his or her 
favorite P{w) and then apply (|TU|) to ( P^D to obtain p{x,y). 

We can consider the effect of dilution. Let < . . . >d denotes averaging with respect 
to Pdiw) defined in (|T^). Then {v^)d = i^ — d) {v^)- Thus T,d{z) is obtained to third 
order in v from the T,{z) in (|3l|) by replacing the moments (v^) with {v^)d- We can 
thus in principle obtain p{x, y) for any value of the dilution. In particular, for — > 1, 
we obtain the simple but uninteresting result 



G{z) = Go{z - (1 - c?) {v{z))) (36) 

The ellipse in (|T^ is distorted slightly. 

For the Cauchy distribution, things simplify enormously. Recall that if the random 
variable x obeys a Cauchy distribution, then (x + constant) and x~^ also obey Cauchy 
distributions. Thus Vk obeys a Cauchy distribution and {v^) = {vY for any integer 
power p. (Note that in contrast to the moments {v'^) the moments {w^) do not exist. 
See (^2]).) With {v"^) = {v)"^, the expression in ( pT] ) suggests that J2{z) is given to all 



w 



1 — wGq 



(32) 
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orders simply by = ' ^^^^ '^'^^ sketch a proof that this is indeed 

the case. 

Indeed with the Cauchy distribution the hermitean version {h 0) of our prob- 
lem has long been known as the exactly solvable Lloyd model [|lT|. The key to the 
solvability is the fact that the averaging integrals over the w^s can all be done using 
Cauchy's theorem. For z in the upper half plane, the integration contours in w can 
be closed in the lower half plane and w is effectively set equal to —27, and vice versa 
if z is in the lower half plane. We must now show that this procedure, proved in the 
literature for hermitean Hamiltonians, continues to hold for non-hermitean Hamilto- 
nians. Imagine expanding (^) in a series in the v^s. Upon averaging, the v^s are 
effectively replaced by {v). Using 

- 1 - (37) 



1 + {v) GoPk 1 + {v) Go 

we obtain from (|29|) that indeed 



E(^) = T-7T^ (38) 
1 + {v} Go 

We see from ( P^D that is equal to =F^7 in the upper and lower half plane respec- 
tively. We thus obtain from ( pO]) 



G{z) = Goiz + i-f) e (Im z) + Goiz- i^) 6 (-Im z) (39) 

The density of eigenvalues again consists of two arcs and two wings. When ^ in 
([10|) acts on Go in ( p9|) we obtain 



p{x,y) = po{x,y + 1) 9{y) + po{x,y --i) e{-y) (40) 

In other words, the two arcs of the ellipse are pushed towards the real axis by a 
distance 7. When in (|TDp acts on the step functions in (^) we obtain the density 
on the two wings 
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V27r B[x,'y,t) 
where we have defined for convenience 



B{x,-f,t) = ^(x2+72)2 + 2t2(72_a;2) + ^4 (42) 

and 



a;mm(7) = (Y(tsinh/i)2 - 72) tanh /i (43) 
The critical value of 7 at which the arcs disappear is given by 

7c = t sinh h (44) 

All states are now localized. Evidently, the geometrical construction described after 
still works. In the one way limit, all these quantities reduce appropriately and 



p{x, y) tends to (|0D and (|2T|) 



In this model, 7ci is again O"*". The wings disappear as 7 ^ 0. We see that for 
insufficient non-hermiticity, h < h^. where 

K = sinh-^ (^) , (45) 

all states are localized. Note that in accordance with Anderson et al for any 
finite non-zero 7, no matter how small, there are localized states for h small enough. 
Trivially, we can express (^4[) in yet another way, by saying that for a given non- 
hermiticity h and site randomness 7, we need the hopping to exceed a critical strength 



^ (46) 



sinh h 

before states become delocalized. 
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3 Bounds in the One Way Limit 



In this section we revisit the one way hmit and extend the analysis in FZ. In the one 
way hmit, the determinant of z — H simphfies drasticaUy to (we take N to be even 
for definiteness) 

det {z-H)= (j[{z - w,)^ - 1 (47) 

(For arbitrary t and h the corresponding formula is considerably more complicated.) 
Note that ( ^7|) is completely symmetric in the {wi}, and thus, for a given set of 
site energies it is independent of the way the impurities are arranged along the ring. 
Differentiating the logarithim of det (2; — H) with respect to z we obtain formally 




Defining 



9{^)^{(^)) (49) 
\z — wj 



and 



1 /d^'-' 



we obtain 



Giz) = J:9i+ii^)i9i{z)r-' (51) 

1=0 

We are to evaluate this in the large N limit and then compute the density of eigenval- 
ues by p{x,y) = ^-^G{z) For an arbitrary P{w) it is non-trivial to evaluate gi {z) 
and to carry out the sum in (^l]). 



In FZ it was noted that in the case of the Cauchy distribution P{w) this program 
can be carried out. It is worthwhile, so that we can see how the expression in (|5l|) 
works, to use (|51| ) to give a somewhat more compact derivation of (pO| ) and (|2T| ) here 
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than in FZ. The key is that in the Cauchy case gi{z) is equal to g{zy, a simple power 
of 



1 



(52) 



z + i'y sign [Im z] 



a result which can be obtained easily by using Cauchy's formula for contour integra- 
tion. The sum in (^) is now immediate, giving G{z) = iJg^^^N (where we neglect the 
difference between and — 1). Thus, for | g{z) \< 1, we have G{z) = g{z), and 
for \ g{z) \> 1, we have G{z) = 0. We obtain 



Differentiating according to ([T0|) we obtained the density of eigenvalues in (|20D and 



While we are not able to evaluate the sum in (0) for P{w) other than the Cauchy 
distribution, we can devise bounds on the domain over which the density of eigenvalues 
is non-zero. Indeed, the analysis given in the Cauchy case suggests how to proceed. 
We will assume that P{w) has bounded support; an example is the box distribution. 
We will always take P{w) to be even for simplicity. Note that gi{z) = gi{—z) then 
and p{x,y) = p{—x,—y) = p{x,—y) = p{—x,y). This allows us to restrict the 
subsequent analysis to the first quadrant {x > and y > 0). 

We begin by noting that if | gi* \< 1, then the /* term in the sum in ( |5TD can 
be dropped. This inequality defines a domain D {x,y] I*, {r}) in the complex plane. 
(Here {r} denotes the set of parameters in P{w).) The intersection of all such do- 
mains D {x,y;l*,{r}), as /* ranges over all integers from /* = 1 to infinity, defines 
a domain D over which the density of eigenvalues p {x, y) vanishes and which may 
of course be considerably smaller than the actual domain over which p {x, y) van- 
ishes. Note the I = term in (|5T|) can of course never be dropped. Ostensibly, it 



G{z) 



g{z)9 {1- \ g{z) \) 




(53) 



(|2ll) with its arcs and wings. 
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contributes to G{z) a term g{z) and hence, according to (0), to p{x,y) the term 
{6{x — w) S{y)) = P{x) S{y). For P{w) of bounded support, this represents a finite 
fine segment and as long as it does not protrude into D it is not relevant to our search 
for a bound on D. Whether or not such a line segment actually appears in p{x,y) is 
an issue upon which we will remark later. 

Given the preceding observation, it is easy to obtain some bounds. Our analysis 
to follow will be rather informal; however, we believe that all the bounds we obtain 
could be made mathematically rigorous. 

As a simple first example of a bound, consider 

I gi{iy) I < IdwP{w) I \' 

which is less than 1 for all positive integers / if | y \> 1. Thus, we conclude that 
for I ?/ |> 1 the density of eigenvalues vanishes. This indeed holds true for all the 
examples worked out in FZ. 

As a second example, we note that 



gi iz)\ = / dwPiw) I —. r (55) 

X + ty — w 



dwPiw)\l- (56) 

\[x — wy + y^ ) 

By definition, x — w > x — Wmax, and hence {x — wY > (x — Wmax)"^ if a; > Wmax- 
We obtain 

\9iiz)\ < (t — ^1 (57) 

\{X- WmaxY + r/ 

and thus there is no density of eigenvalues if 

{X - Wmaxf + > 1 (58) 

and 
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max 



(59) 



Combining this with the result obtained above, we conclude that the density of 
eigenvalues is confined to an oval-shaped region defined by y > 1 and (|58| ) and (|59D . 
(In drawing a figure, it is useful to consider the cases Wmax > 1 and Wmax < 1 
separately.) In particular, we conclude that the tips of the wings cannot extend 
farther than {wmax + !)• (The Cauchy distribution does not have bounded support 
and indeed the wings extend to infinity.) 

The bound (|55D can be improved if P{w) vanishes for —Wmin < w < Wmin, in 
other words, if there is a hole in the middle of the distribution. The sign model 
defined by (|16]) provides an example of this class. In the integral defining gi{z), we 
have (x — w)'^ > (x + Wminf' for w < and (x — w)^ > (x — Wminf for w > 0. We 
find that the condition \ gi{z) \< 1 leads to 



In particular, at x = we find that p{x,y) = when > 1 — Thus, for 

> 1 the entire y axis is free from eigenvalues, and the eigenvalue distribution has 
split into, by symmetry, at least two separate "blobs." This was found to be indeed 
the case in the sign model in FZ with the critical value = 1. (See figure (2c) in 




(60) 




2s' > 1 + r 

(s^ 4:WjYiin-^) ^ 

The inequality is satisfied for (s^ — Awminx) > 1 and thus for 



(61) 




(62) 



There is a similar inequality for x < obtained by reflection. 
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4 Sign Model and Expulsion of Eigenvalues 



We now extend the analysis given in FZ of the sign model in the one way limit. For 
this model, using (^) we obtain 

(tr log - if)) = 2-^^ log [(2-r)"(2 + r)^-"-l] . (63) 

n=0 / 

This expression has the simple physical interpretation that for a specific realization 
of H, n of the impurity potentials have the value +r while the others have the value 
— r. Define n = ^ (1 + '^)- Then for this specific realization, the eigenvalues are 
given by 

(^2 _ (i^)- = e^' (64) 

Here j = 1, . . . , A^, and the N —>■ oo limit is indicated with 6 running between and 
2tt. Thus, the density of eigenvalues p{x,y) spans a one dimensional curve for each 
realization. For finite A^, the density of eigenvalues is formed by superposing these 
curves (See figure (2) in FZ.) In the limit — > oo, the binomial coefficient appearing 
in (p3D is sharply peaked around n ~ A^/2, and can be approximated by e ~ . Thus, 
the range in cr vanishes like (This is consistent with the width of the apparently 
two-dimensional regions obtained by diagonalizing finite A^ Hamiltonians in figure (2) 
in FZ.) 

Thus, we see that the original unit circle spectrum of the deterministic Hq is 
distorted by randomness into the curve 

^2 = ^2 ^ ^ < < 27r (65) 
in the complex z plane. Clearly, Tc = 1 is a critical value of r. For r < 1, the curve 



(|65|) is connected, enclosing a region free of energy eigenvalues. For r > 1 it breaks 
into two disjoint symmetric lobes that are located to the right and to the left of the 
imaginary axis. We mention in passing that if we substitute z —>■ iy in (0) we obtain 
the solution 
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= l-r^ (66) 

independent of a even for finite A^. This is consistent with figure (2) in FZ. 

We mention that our resuh for the sign model satisfies the bounds we just de- 
rived, of course. For instance, the curve (|65|) intersects the positive real axis at 



Emax = A/rM-T (and at Emm = Vr^—l if r > 1.) The bound that the tip of the 



wing cannot extend beyond {wmax + 1) is certainly satisfied since (r + 1) > vrM-T 
for r positive. Note also that by our discussion there are no wings in the sign model. 
This can also be argued heuristically from (|51| ) since term by term there is no cut in 
gi{z). Thus, the critical strength of randomness at which localized states first appear 
(called 7ci in the Introduction) is effectively infinite. 

One feature exemplified by our result for the Cauchy model and the sign model is 
the presence of a "hole" in the spectrum, that is, there are no eigenvalues near z ~ 0. 
In the one way limit we can understand this feature analytically. Rewrite (|^) as 



— tr G{z) 



E n 



Z - Wj, 



(67) 



iV ^ ' V^{z-w,)-\ 

i 

The expression in the square bracket is a polynomial. Thus, near z ~ 0, Gi^z^ is 
singular only if Hi'^^i = 1- So the presence of the hole is related to the probability 
that rii'U^i = 1- This provides another understanding of the critical transition at 
Tc = 1 in the sign model: the probability just referred to becomes \ at that point. 
Our simple argument suggests that the presence of the hole is generic for any P{w) 
which does not have any singularity at w = 1 (such as the box distribution). See 
figure 1 of FZ. 

Earlier, we mentioned that the I = term in (^Tj) apparently contributes a term 
P{x) 6{y) to p{x, y), representing a finite line segment for P(w)'s with bounded sup- 
port, such as in the box model, or representing two points, such as in the sign model. 
From the argument just given and from the explicit result ( |65D this contribution does 
not seem to be there. We believe that this expulsion of eigenvalues from the region 
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around z ~ is a generic phenomenon. 

In order to have a clear understanding of this phenomenon, it is perhaps better to 
give a simple explicit demonstration, rather thatn a heavy formal proof. We consider 
the simplest possible situation, combining the one way limit and the single impurity 
limit. Schrodinger's equation reads 



{E -w)iJi = ^2 (68) 
E^P2 = (69) 

EijN-i = i^N (70) 
E^N = ^1 (71) 

In the absence of the impurity, the spectrum is clearly given by the unit circle | |= 1. 

Suppose in a particular realization of the random energy w we have w > 1. We 
see that there is a solution tp^ = ^ V'l ; i^N-i = ^ V^i , • • • , which decreases 
rapidly, so that by the time we reach the site 2, ?/'2 is exponentially small. Thus, we 
have a real eigenvalue E = w to exponential accuracy. A similar discussion can be 
given for w < —1. In contrast, in a realization in which | w |< 1, the solution just 
described becomes unbounded and must be excluded. 

As we approach the circle moving along either of the two wings (which have 
spectral weight -^), in other words, for E real and tending towards ±1, the localization 
length diverges as 

LiE) ~ , (72) 
log \ E \ 

It is tempting to conjecture that in general, if Ec is the point at which the (right) 
wing attaches to the complex spectrum (for example, Xminil) in (^31) ), the localization 
length diverges as 
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L{E) ~ ^ (73) 

as E ^ Ec- 

Finally, we mention for future work that given (|48|) we can write down easily in 
the one way limit a formal expression for the correlation function 

1 111 



much studied in recent years in the context of random matrix theory. Define 



z — w z — w 



9iz) -g{z') 



z — z 



(75) 



and 



(Of course can be obtained from h by differentiation.) We then obtain 



G(^,^') = Trim Vi/+i(^'^')(V(^>2^')) 



'\\N-1 



+ {N- l)h,^,i> {z, z')hi,i,+,{z, z') {hiv{z, z')f-^] (77) 
For an arbitrary P{w) this is certainly no easier to evaluate in closed form than it is 



to evaluate G{z). (It is however simple to show that in the Cauchy case G{z, z') has 
no connected piece for z and z' on the same side of the real axis.) 
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5 Conclusion 



We have studied the spectrum of a non-hermitean hopping Hamiltonian, with real 
eigenvalues corresponding to localized states and complex eigenvalues corresponding 
to extended states. Going beyond an earlier work FZ, we treat multiple scattering, 
obtaining a general formula, which we evaluated to a certain finite order, arriving 
at an expression ( |35D from which one can compute the density of eigenvalues for an 
arbitrary probability distribution P{w) of the random site energies. For the Cauchy 
distribution, we obtain explicit analytic results. In the so-called one way limit, while 
it is not possible to give explicit analytic results for an arbitrary P{w), we are able 
to obtain bounds on the spectrum of the Hamiltonian. The spectrum of the so-called 
sign model was worked out explicitly and shown to satisfy the bounds derived. We 
also give an analytic argument for the expulsion of eigenvalues from the region around 
z ~ 0, a phenomenon we believe to be generic. 

Note: While completing this paper, we noticed the paper of Goldsheid and Kho- 
ruzhenko |]10[ who also treated the Hamiltonian in (|I|). 
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